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The dispersive interaction between nanotubes is investigated through ab initio theory calcu- 
lations and in an analytical approximation. A van der Waals density functional (vdW-DF) 
[Phys. Rev. Lett. 92, 246401 (2004)] is used to determine and compare the binding of a pair of 
nanotubes as well as in a nanotube crystal. To analyze the interaction and determine the impor- 
tance of morphology, we furthermore compare results of our ab initio calculations with a simple 
analytical result that we obtain for a pair of well-separated nanotubes. In contrast to traditional 
density functional theory calculations, the vdW-DF study predicts an intertube vdW bonding with 
a strength that is consistent with recent observations for the interlayer binding in graphitics. It 
also produce a nanotube wall-to-wall separation which is in very good agreement with experiments. 
Moreover, we find that the vdW-DF result for the nanotube-crystal binding energy can be ap- 
proximated by a sum of nanotube-pair interactions when these are calculated in vdW-DF. This 
observation suggests a framework for an efficient implementation of quantum-physical modeling of 
the CNT bundling in more general nanotube bundles, including nanotube yarn and rope structures. 



PACS numbers: 61.50.Lt, 61.46.-w,71.15.Mb 

I. INTRODUCTION 

Carbon nanotubes (CNTs) have a wealth of exciting 
physical properties that have made them the focus for 
a very broad range of fundamental-science studies W The 
CNTs have, for example, an exceptionally large Youngs 
modulus P The individual CNTs have nanoscale diam- 
eters and micronscale lengths but a range of CNT as- 
sembly processes promise technology applications even 
on more macroscopic scales. Thermal treatment can 
cause a fullerene source to transform into a highly regular 
CNT crystal with parallel tubes aligned in a hex agonal 
structure.^ The tubes can also form CNT bundlesS^E 
in which essentially parallel CNTs still have a very high 
deg ree of (local) order. The bundles can be spun into 
yariJS12l and further twisted to produce torque-free ropes 
of micrometer diameter (and arbitrary length). The yarn 
and ropes have a large strength and a unique ability to 
absorb elastic energy in reversible extensions.^ By pre- 
selecting the nanotube source material, for example, as 
single-walled CNTsj^ it is possible to ensure specific 
physical properties (like m etalli c conductivity) also of the 
resulting well-aligned yarnPSl 

The sciencePEHIIlE] and technology progressing!!] 
challenges us to present a quantum-physical characteri- 
zation of the bonding in the nanotube crystaP and, by 
extension, in the bundles. It is valuable to have a method 
for parameter-free characterization of general CNT bun- 
dles and it is important to test the accuracy of avail- 
able computational tools. The CNT crystals (and bun- 
dles) are (approximately) periodic and have a relatively 
simple order. This makes them accessible to calcula- 
tions in density functional theory (DFT) which, in princi- 
ple, provides quantum-physical accounts of general mate- 



rial bonding. It is straightforward to provide quantum- 
physical calculations for a parameter-free characteriza- 
tion of the intra-CNT electronic and atomic organiza- 
tion using traditional implementations of state-of-the-art 
DFT calculations! 12 * 13 ! These calculations use either the 
local density approximation (LDA)M^ and/or the semilo- 
cal generalized gradient approximation (GGA), for ex- 
ample, as parameterized in the Perdew-Burke-Ernzerhof 
flavor. 15 However, the CNT crystals and bundles are 
graphitic materials and the intertube attraction is known 
to be dominated by relati vely s oft dispersive or van der 
Waals (vdW) interactions.^ Neither LDA nor GGA 
provide any ph ysics-based a ccount of the bonding be- 
tween the CNTs! 17 ' 18 ' 19 ' 2 ^ 

In this paper we use a rec ently developed van der 
Waals density functional 22 ^! (vdW-DF) to provide a 
quantum physical account of the vdW bonding in a 
hexagonal crystal 3 of parallel semiconducting nanotubes. 
We perform state-of-the-art DFT calculations of the 
intra-nanotube structure within GGA and of the inter- 
nanotube binding within vdW-DF. The study testifies to 
the strength of this vdW bonding, which is normally de- 
scribed as soft but nevertheless contributes significantly 
to the cohesion of the CNT crystal. Our results allow 
a test of the vdW-DF theory method by comparison 
against structure measurements for the highly ordered 
CNT crystaP and bundle^ structures. The study supple- 
ments a recent vdW-DF calculation^ on a simple poly- 
mer, polyethylen, for which there also exists experimen- 
tal characterization of the crystaline structure.^ It also 
supplements vdW-DF calculations of the benzene and 
DNA base-pair interactions 2 ^ in a wider program on cal- 
culating dispersive interactions in carbon and organics 
materials. We provide a parameter- free theory determi- 
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nation of the CNT bonding in the crystal and between 
pairs of parallel nanotubes and document that a summa- 
tion of nanotube-pair interaction energies (calculated in 
vdW-DF) represents a fair approximation for the vdW- 
DF results for the crystal bundling energy. We also detail 
the nature of the mutual CNT interactions by identify- 
ing a set of distinct vdW interaction regimes. We show 
that the vdW interaction is significantly enhanced at the 
bonding separation compared with the value estimated 
from the asymptotic interaction. 

The vdW-DF calculations correct the accuracy prob- 
lems arising in traditional state-of-the-art implementions 
of DFT (that uses LDA or GGA) without loss of the 
traditional-DFT scaling 23 (computation cost increasing 
~ 0(N 3 ) with system size). DFT calculations in GGA 
show no meaningful binding. 27 While LDA calculations 
can mimic the CNT binding it underestimates the bind- 
ing separation. For the CNT bundles, the LDA resultP^l 
for the wall-to-wall separation, Ale>a = 3.1 A, is 10% 
shorter than the experimental value, Acnt = 3.4 A. 
Moreover, the LDA result,^ « 10 meV/atom, for the in- 
tertube binding in a crystal of metallic (6,6) CNTs is sig- 
nificantly smaller than the estimate^ (50 meV/atom for 
graphitics materials) extracted from measurements of the 
binding of polyaromatic hydrocarbon (PAH) molecules 
on graphite. The vdW-DF method corrects those prob- 
lems without loss of scaling advantages by supplementing 
the LDA for correlation with a nonlocal contribution^ 
that scales like 0{N 2 ) with the system size.^The vdW- 
DF method clearly has a better scaling than implemen- 
tations of canonical Moller-Plesset perturbation theory 
(MP2) for extended structures like polymer crystals.^ 
Specially adapted MP2 implementations can achieve a 
linear scaling with size for large molecules.^ The adapted 
MP2 method has also been applied to extended one- 
dimensional systems.^ It is unclear how the adapted- 
MP2 evaluation of correlation and the vdW-DF determi- 
nation of nonlocal correlation compare in actual comput- 
ing cost for extended structures like polymer cr ystaliP 4 ^* 
and for large bulk and surface-adhesion systems! 21 * 33 * 34 * 35 -^ 
In any case, the complete MP2 calculation^ also involves 
Hartree-Fock calculations which scale worse than gen- 
eral DFT implementations (including, for example, vdW- 
DF). 

The paper is organized as follows. In section II we 
identify regimes of interactions in crystals of nanotubes 
and discuss qualitative differences in the vdW bonding 
of semiconducting and of metallic nanotubes. Section III 
presents a summary of the vdW-DF calculation method 
that we use to obtain an ab initio characterization of 
the semiconducting-nanotube crystal. In section IV we 
present an analytical description of the nanotubc inter- 
action and in section V we discuss both the strength and 
nature of the inter-nanotube interaction. Section VI con- 
tains our conclusions and acknowledgements. 




FIG. 1: Atomic and electronic structure as well as filament 
organization in crystals of semiconducting (8,0) zigzag nano- 
tubes. The lower panel shows the fully relaxed atomic con- 
figuration of the individual (8,0) nanotubes as calculated in a 
traditional implementation of ab initio density functional the- 
ory (DFT). The top left panel shows our corresponding tra- 
ditional DFT results for the length-averaged electron density 
concentrated at the atomic nuclei; the contour spacing is spec- 
ified in steps of 0.15 eA 3 . Finally, the top right panel shows 
the hexagonal crystalline order of the nanotube bundle. We 
calculate the intertube dispersive interaction and determine 
the nanotube crystaline structure using a recently developed 
ab initio van der Waals density functional approach!^! 



II. REGIMES OF NANOTUBE INTERACTIONS 

There are several regimes of interactions relevant for 
the cohesion of the nanotube crystals. The individual 
tubes are held together by exceptionally strong covalcnt 
bonds between neighboring carbon atoms separated by 
just 1.42 A. The binding between the CNTs is instead 
dominated by the vdW interaction that binds the tubes 
at a wall-to-wall CNT separation of 3.4 A. The vdW 
interaction also causes an inter-tube attraction even at 
asymptotic distances. There are qualitative differences 
in the vdW forces in the asymptotic regime (where the 
interaction is defined by the dipolar electrodynamical re- 
sponse) and at bonding separation in graphitics materials 
(where multipolar contributions arc documented to en- 
hance the interaction. )p2 

Moreover, the vdW interaction between extended 
semiconducting and metallic structures (for exam- 
ple, semiconducting and metallic CNTs) is qualita- 
tively differ ent, at leas t in the regime of asymptotic 
interactions ! 23 | 36 | 37 | 38 | 39 ] Single- walled CNTs can exhibit 
both a semiconducting and metallic nature of conduc- 
tion depending on their chirality. Being low-dimensional 
systems, the semiconducting and metallic CNTs there- 
fore exhibit significant differences in their electronic 
response a nd, conseq uently, in the asymptotic vdW 
interactions! 36 * 37 * 38 * 39 * In the strictly asymptotic regime it 
is possible to view the CNT as wires. For a pair of insulat- 
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ing or semiconducting wires, the asymptotic form of the 
interaction is known to eventually aquire a d~ 5 depen- 
dence with the separation d between the wire centers of 
mass. In contrast, Dobson et al. recently used a coupled- 
plasmon expansion and approximations valid for asymp- 
totic mutual separations to derive a mutual interaction 
energy with a —d~ 2 (log(d)) -3 / 2 asymptotic scaling for 
metallic wires It is not known to what extent qualita- 
tive differences between the vdW binding of metallic and 
semiconducting nanotubes persist down to distances rele- 
vant for their binding in CNT pairs or bundles^ but that 
is an important question beyond the present scope. Ex- 
traction of CNT binding energies from the metallic-wire 
study Ref. [39] is complicated because there are two con- 
voluted interaction effects arising as the CNTs approach 
each other. Firstly, the C NT morphology (a hollow cylin- 
der) manifests itsel f 40 * 4 ^ even when the long-wavelength 
form of the mutual dielectric response remains applica- 
ble. Secondly, the nature of the van der Waals interaction 
changes^ so that it is no longer dominated by the long- 
wavelength response form, but also retains interaction 
contributions defined by the multipole response 

Figure [l] shows schematics of the electronic, (intra- 
tube) atomic, as well as (inter-tube) crystalline ordering 
(bundling) of nanotubes. We study the mutual binding of 
pairs and bundles of CNTs that have a chirality vector-^ 
(8,0). These CNTs have a diameter a little larger than 0.6 
nm, a four-fold rotational symmetry, and an along-axis 
structure repeating itself every 32 atoms. Confirming 
also previous investigations,^ we find that state-of-the- 
art DFT calculations using GGA provide an excellent ac- 
count of the intra-nanotube structural organization and 
electronic properties such as the nature of conduction. 
We use the traditional DFT-GGA results, obtained in a 
plane-wave implementation,* 4 ^! as the starting poin t for 
vdW-DF calculations of the intcrtube binding. 22 23 

Recent d ensity-functional approxim ations by 
ugMMMMM and b y therP™aillSI extend tra- 
ditional DFT to provide a seamless, parameter- free 
characterization of the vdW binding without introduc- 
ing double counting at separations with finite ov erlap 
of electron densities. In our vdW-DF method 2 ^ 2 ^ 
we extract the exchange from GGA calculations but 
supplement the local density approximation for the 
correlation energy by a nonlocal correlation energy 
contribution Efr. This contribution is evaluated from 
the electron densities of the underlying traditional 
DFT calculations in GGA. This vdW-DF description 
remains applicable and effective even for large extended 
systems that are accessible for standard ab initio DFT 
calculation (although at an increase in computing cost). 
In fact, the vdW-DF method exhibits the same scaling 
as the underlying traditional DFT calculations!^ 3 -* The 
vdW-DF approach permits ab initio characterizations of 
large bulk systems, for example, produced by potassium 
intercalation!^ It further permits ab initio investigations 
of very large surface-adhesion systemSj for example, 
graphite-adsorption of PAH molecules ! 33 * 34 * In a con- 



trolled approximation it even permits an ab initio study 
(in a repeated unit cell containing 146 atoms) of the 
adhesion of graphite sheets on SiC surfaces.^ 

The vdW-DF evaluation of the nonlocal correlations 
-E" 1 (vdW interaction energy) involves a density-weighted 
integration of a kerneP^l that co ntains a r ich account 
of the complex electrodynamics.^ * 22 * 23 * 44 * Our vdW- 
DF is not developed to include an explicit account of 
the asymptotic interaction between extended metallic 
one-d imens ional systems.^ The form of the vdW-DF 
kerneP^EHl ensures the correct asymptotic behavior of 
vdW interactions for atoms, molecules, and most surface 
and bulk systems. It also describes the asymptotic in- 
teractions for extended low-dimensional systems that are 
isolating or semiconducting. The form of the vdW-DF 
kernel ensures the correct asymtotic form of the interac- 
tion between pairs and within crystals of the (8,0) CNTs 
because these are robustly semiconducting (characterized 
by a significant gap^). More importantly, our vdW-DF 
calculations of the (8,0)-CNT binding not only eventually 
reproduces a d~ 5 dependence in the interaction energies 
but reveals a much finer structure and remains fully ap- 
plicable at general separations. 

To interpret this rich structure in CNT binding-energy 
variation with CNT separation, we also present in this 
paper an analytical evaluation of the CNT-pair interac- 
tion, valid at intermediate-to-asymptotic distances. Our 
analysis tool (but not our full vdW-DF calculations) 
makes assumptions of non-overlapping electron densities 
and of a long- wavelength form of the CNT dielectr ic re- 
sponse but it respects the CNT morphology.* 41 * 50 *^ Com- 
parison with the full vdW-DF calculations therefore al- 
lows us to deconvolute effects arising from the change in 
nature of the dielectric response. We thereby provide 
an analysis that splits the interaction into two major 
regimes: (1) a close regime at or near binding separa- 
tions where full ab initio vdW-DF calculations are essen- 
tial for an accurate account and (2) an intermediate-to- 
asymptotic regime where the long-wavelength dielectric 
response remains applicable but where the CNT mor- 
phology specifies the variation of the interaction. The 
result of this analys documents that the vdW interac- 
tion enhances at bonding separations compared with es- 
timates that can be extracted from knowledge of the 
asymptotic form of the inter-nanotube bonding. 



III. COMPUTATIONAL METHODS 

The top right panel of Fig. [T] shows a schematic of the 
repeated two-dimensional hexagonal array of the nano- 
tube bundles. We apply the vdW-DF^ method to in- 
clude the dispersive interaction within the framework of 
traditional plane-wave DFT both for dimcrs of nanotubes 
as well as for an infinitely extended nanotube crystal 
(Figure [TJ. A self-consistent formulation of vdW-DF has 
recently been derived, implemented, and tested.^ Here 
we use the original, non-selfconsistent (post-GGA) imple- 



4 





R geo 


lit 


R 




L 



FIG. 2: Length- and radially-averaged electron density n(s) 
shown together with the positions of the effective and the 
geometric radius. The radial separation is given relative to 
the geometric radius of the nanotube. The background inset 
shows the length averaged electron density with contour lines 
separated in steps of 0.15 eA 3 . 



mentation^! that rests upon and utilizes traditional DFT 
calculations to obtain the electron density variations. 

In the present vdW-DF study we furthermore take ad- 
vantage of the success of the traditional semi-local (GGA) 
density functionals to describe the infra-molecular prop- 
erties of the nanotubes (as well as the electron densi- 
ties). It is, in principle, possible to provide an all-vdW- 
DF characterization of the intra-nanotube atomic struc- 
ture (allowing relaxation under vdW-DF forces^) but 
the computation costs would be large. Our previous ex- 
perience from an all-vdW-DF characterization of a sin- 
gle graphite sheetPS suggests that only minute differences 
would result for the CNT structure if we replaced the 
GGA intra-tube characterization by a full vdW-DF char- 
acterization. 



A. Nanotube atomic and electronic structure 

In the actual vdW-DF implementation, a large set of 
state-of-the-art traditional DFT calculations determine 
the atomic structure (of the individual nanotube) and 
the electron density variation of the nanotubes when iso- 
lated, when assembled into a hexagonal crystaline struc- 
ture, or when aligned as a parallel nanotube dimer. We 
use a plane- wave code^l with ultra-soft pseudopotentials 
and a 1 x 1 x 8 Monkhorst-Pack — fc-point sampling of 
the Brillouin zone of the periodically repeated unit cell 
containing 32 atoms per nanotube. We perform self- 
consistent calculations in the Pcrdew-Burke-Ernzerhof 
(PBE)-fiavop^of GGA for the exchange and correlation 
functional. We choose a plane- wave energy cut-off 476 
eV and specify the Fast Fourier Transform (FFT) grid so 
that the (density) grid spacing remains smaller than 0.14 



A in all directions. 

We first determine the atomic structure of the indi- 
vidual (8,0) nanotube, lower panel of Fig. [I] We use our 
ab initio calculations of the strong intra-nanotube atomic 
forces to relax the morphology to a total residual atomic 
force 0.05 eV/A per unit cell. The structure is character- 
ized by an effective geometrical radius i? gC o specified as 
the average distance of the carbon nuclei from the axis 
defined by the nanotube center of mass. We determine 
the value of that geometric radius to be R gco = 3.18 A. 
As explained above, the atomic structure (obtained in 
PBE-GGA for the isolated nanotube) is kept frozen in all 
subsequent calculations (for nanotube crystal and dimer 
cases). We have explicitly tested that no additional in- 
tratube atomic-relaxation is relevant at (or beyond) the 
separation that characterizes the nanotube-crystal/dimcr 
binding in our subsequent post-GGA vdW-DF character- 
ization. 

Next we determine the electronic structure from tra- 
ditional DFT calculations for the isolated nanotubes, for 
the nanotube bundles, and for the dimer structures. The 
upper left panel of Fig. [I] shows contours (at 0.15 eA -3 
interval) in the isolated-nanotube electron density varia- 
tion averaged along the nanotube axis. The maximum in 
electron density coincides with the radial position of nu- 
clei and naturally identifies the nanotube wall. To char- 
acterize the inter-nanotube vdW-dominated binding we 
further calculate the electron density for nanotubes in 
hexagonal crystals and dimers as a function of the wall- 
to-wall nanotube separation A and as a function of the 
nanotube rotation angle (relative rotation angle in the 
case of the dimer study) . 

We find that there are important electron-density over- 
laps for A < 4 A when the self-consistent GGA elec- 
tron densities differ from a superposition of individual- 
nanotube electron densities. The existence of these inter- 
tube electron-density overlaps has direct consequences for 
the details of the intertube vdW binding. However, the 
electron-density overlap does not reflect the existence of 
any relevant (and physically meaningful) binding arising 
within the GGA calculations themselves (in any of the 
GGA flavors). 



B. van der Waals density functional theory 

Like the LDA and GGA functionals, the vdW-DF is 
defined by approximations for the exchange and corre- 
lation energies. The nonlocal dispersive interactions re- 
sponsible for keeping the bundle of nanotubes together 
are included as a significant extension (correction) of the 
correlations energy in the underlying GGA calculations. 
Specifically, we completely replace the GGA description 
of correlation but use the self-consistent GGA result for 
the electron density to evaluate a new correlation energy 
that includes the nonlocal nature of the vdW binding. 

Our vdW-DF splits the correlation up into local and 
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nonlocal contributions 122143 



LDA 



ET 



(1) 



with the local part approximated i n the LDA. The non- 
local correlation energy is expressed 22 * 43 ^ 



E. 



nl 



du 
2^ 



tr[ln(l - Vx) - ln(e)] 



(2) 



where u is the imaginary frequency, V is the inter- 
electron Coulomb interaction potential, x is the local- 
field density response. The isotropic dielectric function, 
e = tr(l + 47ra))/3, is also specified by the local-field 
density response, x = V • ot ■ V. The nonlocal correlation 
energy is further approximated 



E. 



nl r 



- J drc?r'n(r)0(r,r')n(r'), 



(3) 



through a kernel <p specifie d by a number of sum rules 22 
and physics results.^ 2 * 23 * 53 * The interaction energy ^ is 
consistently constructed to vanish for a homogeneous sys- 
tem. The kernel </> is specified by a pair of local pa- 
rameters <7o( r ) and qo( r ') that depend on the electron 
density and the density gradient. The kernel can be 
tabulated in advance in terms of an effective separation 
D = {(q + q' )/2]\r — r'| and an asymmetry parameter 
5 = (q — q' )/(qo + q'o)- The value of the q 's are cho- 
sen to reproduce the (plasmon-pole) response of a weakly 
perturbed electron gas, Ref. [22] 

With the evaluation of the nonlocal energy contri- 
bution (from underlying GGA calculations of the elec- 
tron densities) we arrive at a vdW-DF total energy 
calculation.^ 



E 



vdW-DF 



= E° + Ef, 



E° — E GGA — e gga 



E 



LDA 



(4) 
(5) 



The GGA energy term £; GGA is here evaluated in the 
revPBE 54 flavor (based on the self-consistent calculations 
for the electron density that we obtain in the PBE flavor 
of GGA). Effectively this amounts to a small adjustment 
of the exchange contribution, which we do to minimize 
any potential artifici al exchan ge binding in the plane- 
wave formalism usedP222ElMl The new (semilocal) en- 
ergy contribution E° represents a modification of GGA 
that, for example, retains a description of the kinetic- 
energy repulsion as well as, for example, covalenP^ or 



C. van der Waals density functional calculations 



intramolecular density variation causes a very large con- 
tribution to E^r that must be carefully subtracted in 
our ab initio calculations of the inter-nanotube interac- 
tion. Moreover, the evaluation of this intra-nanotube E£ 1 
energy is somewhat sensitive to the relative position of 
atomic positions and FFT grid pointsP^Nevertheless, ro- 
bust and efficient evaluation of the inter-nanotube bind- 
ing is possible by simply ensuring that we calculate and 
subtract E™ 1 contributions for the isolated nanotube us- 
ing FFT grid points that closely match those of the com- 
posite system, as further explained in Ref. 1211 

In practice we evaluate the vdW-DF binding in the 
nanotube crystals and dimers by supplementing every in- 
teracting nanotube system by a suitable reference calcu- 
lation of the isolated nanotube. We determine the bind- 
ing in nanotube bundles by adding (at every nanotube 
separation) a DFT calculation of the electron density for 
a corresponding isolated nanotube located in a cell of 
double size in each of the two perpendicular directions 
and on a FFT grid that retains the absolute grid-point 
spacing. 21 In our calculations of E™ 1 , for the underlying 
GGA calculation we use a FFT-grid spacing that is al- 
ways smaller than 0.14 A in any direction; this choice is 
found sufficient given our work to carefully synchronize 
the FFT gridding when we calculate the electron density 
for the interacting and isolated nanotube system. 

The real-space implementation of i?" 1 is simply applied 
to extended systems as graphite^ and polyethylene^ by 
evaluating E^ 1 for the unit cell electron density as well 
as the nonlocal interaction from its surrounding images. 
To accelerate the vdW-DF characterization we limit the 
evaluation of the multidimensional integral ^ to con- 
tributions from grid points having a density larger than 
1CP 4 a.u. The use of such a density cut-off is strongly 
motivated by the excellent convergence that we have pre- 
viously documented for graphitic systems even when a 
significant ionic bond supplements the binding from non- 
local correlations.^The nonlocal correction from the sur- 
rounding electron density converges rapidly in terms of 
the separation to the unit cell. We have tested that it is 
in general sufficient to only include the interaction from 
the electron density that is less than 12 A away from the 
unit-cell boundaries in the direction along the nanotube 
and 15 A in the directions perpendicular to the nanotube 
axis. Nevertheless, to converge the E™ 1 calculations to a 
sub-meV level and retain a very high relative accuracy 
even in the asymptotic regime, we choose to retain E^ 1 - 
contributions originating from points closer than 24 A 
(and in some cases even closer than 30 A) in the direc- 
tion of the CNT extension. 



The evaluation of the nonlocal correlation requires 
extra care due to a grid sensitivity of the functional 
form ([3]). The vdW binding in the nanotube crystal and 
dimer cases arises almost exclusively from a difference in 
nonlocal-correlation energy E™ 1 for the crystaline/dimer- 
structure and for the isolated nanotubes. However, the 



IV. VAN DER WAALS INTERACTIONS AT 
INTERMEDIATE TO ASYMPTOTIC DISTANCES 

From our ab initio vdW-DF calculations of the asymp- 
totic van der Waals interactions we extract an analytical 
determination of the van der Waals interaction energy 
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-E-vdw per unit length L for a nanotube pair as a func- 
tion of the separation d — 2i? goo + A. The analytical 
result for E v ^w{d) / L rests on the approximation sum- 
marized in Refs. 1401411501 It assumes that the electron 
densities of the two nanotubes do not overlap and con- 
stitute a lowest-order expansion^ of Eq. ^ in terms of 
the external-field susceptibilities, a c g: 



E, 



vdW 



du 
2^ 



itr[a e s,iT 12 a e s t2 T 21 ]. 



(6) 



Here T,-j denotes the dipole interaction tensor, Tjj = 
— ViVj|ri — r j | — 1 . The analysis is possible to car ry out 
for nonisotropic external-field susceptibilities j^EH but for 
an interpretation of our vdW-DF calculations of CNT 
interactions it is sufficient to consider isotropic suscepti- 
bilities a c g. We focus on the interaction regime where 
effects of the CNT morpholog y dom inates the variation 
of E'f w -Evdw with distance! 40 1 50 1 We assume a long- 
wavelength form of a ff so that the resulting analyti- 
cal determination remains valid at such intermediate-to- 
asymptotic interaction distances. 

The physics of the local-field and external-field sus- 
ccptibiliti es define s the parameterization of our vdW-DF 
method! 20 * 22 * 23 * 43 -] The long-wavelength electrodynamical 
response determines the interaction at large distances^ 
and our vdW-DF method describes this response by the 
isotropic effective (external-field) susceptibilit}^ 



'off 



(u;r) 



n(r) 



u 2 + [9<? 2 (r)/(87r)P 



(7) 



We stress that our extraction of this long-wavelength 
form serves only to establish formal connection between 
the full vdW-DF calculations and the analytical approx- 
imation. We also emphasize that neither the effective 
response ([7| nor the full vdW-DF response function^ is 
explicitly designed to accurately reproduce, for example, 
the static dielectric response, in contrast to the functional 
approaches described in Refs. I20|41|44l and l40l Rather, 
the full vdW-DF response functional is constructed ex- 
clusively from an ansatz for the plasmon-pole response, 
conservation rules and many-body calculations^2152l to de- 
scribe the average response. The full vdW-DF descrip- 
tion involves contributions from different frequencies and 
wavelengths and it is the average response, rather than 
the long- wavelength limit that determines the inter- 
actions at binding separations where our vdW-DF ap- 
proach is most needed. 

Our analysis focus on the intermediate-to-asymptotic 
separations furthermore allows us to consider the con- 
tributions to the susceptibility from the electron density 
averaged over the angular and along-tube variations. We 
thus substitute a c g(r;u) — > a e ff(s;w), where s denotes 
the radial distance from the nanotube center. The nano- 
tube interaction per unit length, given in terms of the 



effective response (Fl, becomes!^ 
E vr tw f 00 du 



o 27r Jo 



dsisi / ds 2 s 2 (8) 
Jo 

<5off(si; u)a cS (s 2 ; u) G a ^(sx, s 2 ), 

a,/3=s,e,z 



with the geometry factors 

/■2-7T /■2-7T /"CO 

Ga/3 = / d6x d6 2 / d(z 2 - z x ) (9) 
Jo Jo Jo 

[ T ii '( s iA) zi; s 2, fa, z 2 )] 2 , 

where the dipole interaction tensors are expressed in 
cylindrical coordinates. The result (|8| is easily expanded 
in the inverse center-of-mass separation c?" 1 of the nano- 
tubes, yielding interaction energies of the form 



E. 



vdW 



B5 B 7 
d 5 d 7 



with B 5 and B-? given by 



B, = 



9 



duE<- Q \u) 2 , 



(10) 



(11) 







225 



B7 = = Wj duS(0) («)2 (2) (w). (12) 



Here 



S (l) (u)= / ds2Trsa^(s;u)s i , 



(13) 



is simply the i'th moment of the effective response. We 
calculate the coefficients B 57 directly from the a e g (s; u) 
variation specified by our underlying GGA-DFT calcu- 
lations of the CNT electron density variation. We have 
explicitly tested consistency of this asymptotic evaluation 
and the set of full vdW-DF calculations for A > 16-20 A. 

The relevant external-field electrodynamical response 
(a e s) of the nanotubes is dominated by contributions 
at some radius R > R gco . That is the experience 
gained from describing the electrodynamical response 
and van der Waals interactions of surfaces^ and from 
previous i nvest igation of the van der Waals bonding in 
graphiticsP^^H While the results presented in Refs. 140 14T1 
H0|made the assumption that the response a e g arose ex- 
clusively from the atom w all (at i? goo ), the formal special- 
function evaluatio n 1 41 * 58 ! 



E v dw(d, R) _ B 5 
L ~ d 5 



3 r 2 



F 9 , [ i S 1:1,1 



2 ' 2 ' 2 ' 



4R 2 

IF 



(14) 



is possible as long as we may assume the response a c g 
dominated by contributions at any (single) radius R. The 
interaction result ( 14 1 simply reflects the morphology (in- 
teraction of two hollow cylinders) . The effective vdW-DF 
response a^(s; u) is dominated by contributions at -R goo 
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FIG. 3: Nanotube binding energy E vdw - DF (A) - 
_E vdw-DF (A — > oo), (per unit length) for semiconducting 
(8,0) nanotubes evaluated in vdW-DF as a function of the 
wall-to-wall separation A. The top panel shows vdW-DF re- 
sults for the hexagonal crystal and compares the crystal in- 
teraction energy (thick solid curve) against an estimate (thin 
solid curve) based on a sum of vdW-DF results for the nano- 
tube pair interactions. The bottom panel reports vdW-DF 
calculations of the binding of two parallel nanotubes (dashed 
curves) in three different atomic configurations indicated in 
the insert. The sum of those three pair interactions consti- 
tutes the approximation for the crystal interaction (thin solid 
curve) in the top panel. 



V. RESULTS AND DISCUSSION 
A. van der Waals bonding in a nanotube crystal 

The top panel of Fig. [3] reports our ab initio cal- 
culation (thick solid curve) of the vdW-DF binding 
in a hexagonal crystal of semiconducting (8,0) nano- 
tubes. The vdW-DF result for the binding separation, 
Abind ~ 3.45 A, is in very g ood agreement with ex- 
perimental observations] 4 * 5 * 6 ^ A cxp w 3.4 A, and cor- 
rects the poor structure predicted in traditional LDA 
calculations^ Aj^P^ = 3.1 A. The binding energy for the 
nanotube bundle is large, -E^ind = — ^0 meV/atom cor- 
responding to —0.225 eV/A, consistent with interaction 
strengths that we have previously calculated in vdW-DF 
for the interlayer binding in graphite!^ —50 meV/atom. 
The vdW-DF binding energy is significantly larger than 
the LDA result 10 meV/atom, obtained for a metal- 
lic (6,6) nanotube. 

The figure documents differences between the vdW-DF 
calculation for the full CNT crystal (thick solid curve) 
and corresponding approximations based on CNT-pair 
contributions (thin solid curve). The regular vdW-DF 
calculations yield a CNT-crystal binding energy that is 
larger than the binding energy estimate obtained from 
the summation of pair contributions E^nd' 08 * = —29 
meV/atom corresponding to —0.220 eV/A. We find that 
the vdW-DF energy difference ££££ - E"^ cst is split 
evenly between contributions E^ 1 and E . 

Nevertheless, the vdW-DF results for the pair inter- 
action energies constitute a fair approximation of the 
hexagonal ordering arising in the nanotube bundles. It 
is thus possible to use vdW-DF calculations of the CNT- 
pair interactions at general (parallel) configurations (of 
different relative rotations) to model the cohesion and 
binding in more general nanotube structures such as yarn 
and rope. 



and outside in the CNT density tails. In this paper, we 
choose a value for the effective response radius 

and obtain an analytical approximation 

£ c nl (A) w E vdW (d = 2R gco + A, R cS ). (16) 

that exactly reproduces the asymptotic variation of the 
full vdW-DF calculations up to the second spatial mo- 
ment given by B 5 and By. The relative position of R c g 
and -Rgeo are shown in Fig. [2j 



B. van der Waals bonding in a pair of parallel 
semiconducting nanotubes 

Fig. [3j bottom panel, reports our ab initio calculation 
of the vdW bonding between pairs of parallel (8,0) nano- 
tubes at three configurations 'I', 'IF, and 'III', identified 
in the insert. These are the configurations that are rele- 
vant for the pair-interaction estimate of the CNT hexag- 
onal crystal (thin solid curve in top panel). Even for 
a CNT pair, the nanotube binding is very significant, 
-E'bind w —9.2 me V/ atom, but occurs at slightly differ- 
ent binding separations for different relative nanotube 
rotations. We find that the vdW-DF results for the non- 
local correlation term E™ 1 are almost identical (smaller 
than 1% variation outside binding separations) for the 
three CNT-pair configurations. As is evident in the insert 
(which identifies actual atomic organization investigated 
in our vdW-DF method,) the atomic organization is in 
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better registry for some organization than others. There 
consequently exists some electron-density variation with 
the rotations and our vdW-DF method is sensitive to that 
variation since the semilocal contribution Eq contains a 
description of the kinetic-energy repulsion. 

As an interesting aside, we note that the high symme- 
try of the semiconducting (8,0) nanotube permits us to 
test the grid-sensitivity and consistency of the vdW-DF 
calculations. There must exist a four-fold symmetry in 
atomic positions around the (8,0) nanotube and such an 
approximate symmetry also emerges as a result of the 
initial atomic relaxation that we perform for an individ- 
ual nanotube in traditional DFT. The symmetry implies 
a periodicity tt in the variation of the vdW binding be- 
tween a pair of nanotubcs with the relative rotation angle 
0. However, the imperfect relaxation causes small vari- 
ations in the exact atomic location relative to the grid. 
We find that the vdW-DF calculations are more sensitive 
than the underlying traditional DFT calculations. Nev- 
ertheless, the vdW-DF calculations respect the symmetry 
and produce vdW interaction energies for and + 7r 
relative rotations that are identical even at a sub-meV 
energy scale. 



Approximative microsopic modeling for general 
nanotube-bundle structures 



The comparison between the vdW-DF results for the 
nanotube crystal and for the approximation based on a 
sum of nanotube-pair interactions, Fig. [3] top panel, sug- 
gests a framework for an approximative microscopic mod- 
eling for the binding in more general bundles of (semi- 
conducting) nanotubes. A simple mapping of the binding 
energy for two parallel nanotubes for all combinations of 
independent rotations (relative to the interaction line) 
provides the starting point. Adding such general pair- 
interaction contributions allows vdW-DF calculations to 
account for general vdW bonding in aligned nanotube 
structures, including nanotube yarn and ropes. 

Moreover, the finding of insignificant differences be- 
tween the E£ energy contributions for the three 
nanotube-pair configurations investigated in Fig. [3] sug- 
gests an additional speed up in the modeling. Assuming 
that general, independent nanotube rotations also causes 
insignificant i?" 1 differences, it is sufficient to supplement 
one calculation of i?" 1 (detailed below) with a mapping 
of the general E° variation. This can be obtained at 
a computational cost equal to that of traditional im- 
plementations of DFT. A forthcoming study will pro- 
vide vdW-DF results for the bundling of a broader set 
of semiconducting nanotubes and an explicit test of the 
^vdW-DF anc j £jnl var } a tion with general (independent) 
nanotube-rotation angles to detail the suggested approx- 
imative modeling approach. 



> 

CD 



-150, 




-100 



A [A] 



FIG. 4: Nonlocal correlation energy per unit length for a 
nanotube dimer near binding separations (main panel) and 
in the intermediate to asymptotic regime (insert panel) . The 
thick dashed curves show the results of the full vdW-DF cal- 
culation of the nonlocal correlation energy (vdW interaction) . 
The dotted and dashed-double-dotted curves show the (tradi- 
tional) asymptotic interaction estimates as determined from 
the asymptotics —B^dT 5 and —B$d~ 5 — B-jdT 1 , respectively. 
Finally, the dashed-dotted curve shows the analytical eval- 
uation that approximates the electrodynamical response by 
the long-wavelength form but retains a full description of the 
nanotube morphology. The analytical result also reflects the 
surface-physics insight that the electrodynamical response is 
dominated by contributions outside the radius defined by the 
atomic positions, Fig. 2. 



D. Nature of the vdW bonding at bundle and at 
intermediate separations 



nl 



Fig. [4] compares the full vdW-DF calculation of E' t 
contribution to the CNT-pair interaction, thick dashed 
curves, with vdW-DF based approximations -E v dW; dot- 
ted and dash-dotted curves, near binding separations 
(main panel) and in the intermediate-to-asymptotic 
regime (insert). The contribution E^ 1 is evaluated for 
the configuration T shown in the insert of the lower 
panel of Fig. [3] (but i?" 1 exhibits only insignificant dif- 
ferences between configurations 'I', 'IF and 'III'). All of 
the estimates £vdw are, of course, independent of the 
nanotube rotation by construction. The dashed-dotted 
curves show the analytical CNT-pair interaction estimate 
( fl6| that invokes a long-wavelength form of the electro- 
dynamical response but respects the morphology of the 
interaction problem! 40 * 41 ! The dotted and dashed-double- 
dotted curves show (for d = A + 2i? gco ) traditional inter- 
action estimates, —B 5 /d 5 and —B 5 /d 5 — B 7 /d 7 respec- 
tively. The traditional interaction estimates clearly only 
become applicable in a very remote asymptotic regime 
beginning at A > 16 A. 

The main panel shows that the full vdW-DF cal- 
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culations are necessary around the binding separations 
Abind ~ 3.5 A. Here the interaction is significantly en- 
hanced compared with estimates based on the asymp- 
totic dipolar response. The enhancement is consistent 
with the behavior documented for graphite interactions 
as described in an earlier generation of vdW-DF. 43 The 
enhancement relative to the analytical approximations 
persists even beyond separations (A « 4 A) when there 
no longer exists an overlap of electron densities. It arises 
in part because the complete interaction also contains 
multipole interactions. 20 

However, the contrast between the main panel and the 
insert panel in Fig.[4]also documents a qualitative change 
in nature in the mutual interaction with increasing sepa- 
ration. Gradually there is a transition in an intermediate- 
to-asymptotic regime (shown insert panel) where the mu- 
tual interaction is essentially specified by the morphology 
of the nanotube density variation as summarized in the 
analytical interactions estimate (16 1. 



VI. CONCLUSIONS AND 
ACKNOWLEDGEMENT 

We have presented ab initio calculations of the binding 
in nanotube bundles and in nanotube dimers. Our cal- 
culations rest on a density functional description 2 ^ that 
includes accounts of the dispersive forces. We have, in 
addition, presented an analytical evaluation valid at in- 
termediate to large nanotube separation. 

Our microscopic theory of the CNT binding of semi- 
conducting (8,0) CNTs provides a number of results 
based on the ab initio vdW-DF calculations. The CNT 
study supplements recent microscopic theory studies of 
elementents of the DNA base-pair interaction^ and of 
the polyethylene polymer crystaP^ in a broader goal of 
developing a microscopic theory of self-organization and 



bundling of nanoscale filaments. This vdW-DF study 
finds a nanotube wall-to-wall separation in very good 
agreement with experiments and predicts a vdW bonding 
with a significant strength, consistent with recent mea- 
suments for graphiticsP*] 

Our work furthermore constitutes an analysis that de- 
tails the nature of the mutual CNT interactions by iden- 
tifying a set of distinct interaction regimes. We pro- 
vide an analytical approximation for the CNT pair in- 
teractions at distances when the electron densities are 
nonoverlapping and the dieletric response are dominated 
by the long-wavelength form. Comparing against our ab 
initio vdW-DF calculations (valid at general distances) 
we thereby identify a relatively broad intermediate-to- 
asymptotic regime where the interaction form is primar- 
ily defined by the CNT morphology. 

Finally, this introductory study also suggests a frame- 
work for an efficient implementation of quantum-physical 
modeling of the CNT bundling in more general geome- 
tries, including nanotube yarn and ropes. The vdW-DF 
study documents that a summation of nanotube-pair in- 
teraction energies represents a fair approximation for the 
nanotube-crystal binding energy when the CNT-pair in- 
teraction is calculated in vdW-DF. A simple vdW-DF 
mapping of the nanotube-pair interaction for general (in- 
dependent) CNT rotations relative to the interaction axis 
therefore provides adequate input for describing the vdW 
bonding in general aligned CNT structures. 

We thank B. I. Lundqvist for useful discussions and 
L. Glasser for suggesting and detailing the special- 
function evaluation ( 14 1 from our corresponding parallel- 
nanotube interaction result in Ref. 1501 Support from 
the Swedish Research Council, the Swedish Foundation 
for Strategic Research, the Swedish National Graduate 
School in Materials Science, as well as allocation of com- 
puter time at SNIC (Swedish National Infrastructure for 
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